Limits of the Quantum Monte Carlo method 



Z. Neda and Z. Dezso 
Babe§-Bolyai University, Dept. of Theoretical Physics 
str. Kogalniceanu 1, RO-3400, Cluj-Napoca, Romania 
(Last revised 5 Dec. 1999) 



Abstract 



We consider the one-dimensional quantum-statistical problem of inter- 
acting spin-less particles in an infinite deep potential valley and on a ring. 
Several limits for the applicability of the quantum Monte Carlo methods were 
revealed and discussed. We show the inapplicability of the quantum Monte 
Carlo method for ring-like geometries, realize an unphysical frustration for 
interacting fermions and a minus-sign problem for interacting bosons. 
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I. INTRODUCTION 

Quantum Monte Carlo (QMC) [[[]] methods are known to be powerful tools in studying 
various quantum-statistical problems like quantum-spin models f2j, Hubbard or t-J models 
0, quantum chromodynamical systems [f|] or interacting fermions in real space 0]. It is 
believed that the method has great potentials to handle analytically difficult and complex 
problems. The main idea of the method is to transform a d dimensional quantum-statistical 
problem in d + 1 dimensional classical one, by a procedure closely related to the Path integral 
formalism of quantum mechanics ||. Due to the increased dimensionality the method is 
computer-time consuming. Nowadays, the accessible and powerful modern computers makes 
the method more and more popular. 

As an exercise for the applicability of QMC methods we studied the quantum-statistical 
problem of several short-range interacting quantum particles in one-dimensional space. Dur- 
ing this exercise we found however several technical and conceptual difficulties which limits 
the applicability of the method. The present paper intends to discuss the problem in this 
sense. 
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II. THE METHOD 



In this section we review briefly the method of studying short-range interacting quantum 
particles in real ID spaces by the QMC method. We will neglect the spin variables and 
consider idealized spin-less fermions or bosons. 

The time-independent Schrodinger equation for a single particle with mass m in ID is: 

h 2 

= — Aip + V(x) = Erp. (1) 

8n 2 m 

Discretazing the problem on a lattice with sites of length a, the differential equation (0) can 
be written in the form of L (aL = I with I the length of the considered space) coupled linear 
equations 

h 2 

(-5-2 2X^+1 + ^-1 - 2 • + Vi ■ Vi = E^, (2) 

Burner 

where ipi and Vi denotes the medium values in box % of the lattice for the functions ip{x) 
and V(x) respectively. Introducing the notations 

h 2 

t = c~~2 2 3 

Wi = Vi+ — (4) 

equation @ becomes: 

- t ■ + Vi-i) + Wiipi = EiPi (5) 

We can write now this system of equations in a second quantized form, using as state vectors 
the rii occupation number of the cells: | ni,n 2 , ■■■,rii, ...,til >■ We consider the cf creation, 
di annihilation and hi = cfdi occupation number operators defined for bosons as 



c+ I ni,n 2 , ...,n i} n L >= y/n { + 1 | n Xl n 2 , ...,7^ + 1, ...n L > (6) 
di I rii, n 2 , rii, n L >= ^/r?" | n 1 , n 2 , ....n^ - 1, ...n L > (7) 

and for fermions (n« = {0, 1}) as 

cf I ni,n 2 , ...,rij, ...,n L >= (-l) Sl (l - n<) | ni,n 2 , + 1, > (8) 
Cj I ni, n 2 , n i; n L >= (-l) Sl rii \ n 1 , n 2 , - 1, n L >, (9) 

where s, = Z^C^n^. The wave function of the system becomes 

L 

^ = E^A + |o,o,,....,o>, (10) 

and equations @ can be written as: 

L L 

£H(c+ ^ + cfA+i) + Wicf^} I 0, 0, >= £ £^c+ I 0, 0, > . (11) 

i=l i=l 
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The second quantized form of the Hamiltonian in the discretized space will be: 

H = -tJ2icic* + i + ct,A] + E Win*- (12) 

i i 

Considering more then one non-interactive particles in the same lattice the Hamiltonian is 
unchanged. Interactions between the particles will be introduced via extra terms like 

H = V -^2U^-l), (13) 

i 

or: 

#1 = T4- (14) 



Term (|13|) represent an on-site repulsion, and ([14]) the interaction between particles in neigh- 
boring cells. Depending on the sign of V\ this interaction can be attractive or repulsive. 
The partition function is given by 

Z = Tr(exp (-(3 ■ H)) = < n 1 , n 2 , ...,n L | exp (-(3 ■ H) \ n 1 , n 2 , ...,n L >, (15) 

{«.} 

where the sum is over all the possible combination of the occupation numbers subject to 
the J2i=i n i = N (N the total number of particles) restriction and (3 = l//cgT. 

For a practically successful realization of the MC simulation we have to rewrite the Z 
partition function in a form in which calculation of the P transition probabilities are easy 
when only a few rti occupation numbers have been changed. We will sketch how is possible 
this. 

Denoting by Vi the interaction terms of the forms V hi(hi — 1) and, Vihih i+ i, we separate 
our Hamiltonian in two commuting parts H a and H^. 

H a = (-tc+C 2 - tc+Ci + y + y) + 
+ (-tc^c 4 - tc|c 3 + y + y) + 

... + (-tcl^CL - tcJCL^x + + y ) = 

= # 3 + ••• + #£-!, (16) 



with: 

H t = -tc+c i+1 - tcl +1 c t + y + (17) 

In a similar way Hj, = H 2 + + + Hi. One will observe immediately that the terms 

inside H a and Hb commute, and so: 

< {rii} | exp (-/3H a ) \ {rij} >=< {rii} \ exp (-(3Hi) \ {rij} > ■ 

■ < {rii} | exp {-(3H 3 ) | {n 5 } > ■ < {m} \ H L ^ \ {rij} > . (18) 
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The same equation is valid for H b . 

Ideal it would be to write the partition function as product of exponentials, each of them 
containing one Hi term. In this way the calculation of the transition probabilities when only 
a few occupation numbers are changed would be easier. To achieve this we need to write 
exp [—/3(H a + H b )} as: exp (—(3H a ) ■ exp {—(3H b ). Unfortunately this equation is not valid 
anymore, because H a and H b do not commute. 

However, for two A and B sufficiently small operators one can use the approximation: 

e A. e B = e A + B+OMA,B]~ e A + B_ (1Q) 

Considering now M a large integer, so that and is small enough, we are able to use 
the ([TJ]) approximation in the form proposed by Trotter and Suzuki J7J: 

Z=Y< M I (e*P [~ ■ (Ha + H b )]) M | {m} >^ 

{rii} 



Y < {*} | (exp [-jL . Ha ] ■ exp [-jL • H b ]) M \ {n,} > (20) 

{rii} 



In order to write Z as a product, which in a MC simulation can be considered as transition 
probabilities, we insert between each exponentials in (|20|) a complete set of states: 

I n i,i> n 2,j, ■■■■,n L ,j >< nij,n 2 ,j, ....,n L)j \ . (21) 

Because each of these 2M — 1 sets represent an independent sum over all possible states, 
they will be indexed with a new label, j, different from the spatial one, i. The initial 
trace contribute also to the independent set of states, so the j index can take 2M values. 
Remembering that Hi acts only on the states rii an d n i+ i we can write: 

2M L 

2« e nn^ (22) 

When j and % are both even or odd the Py factors are calculable as 

Pi,j =< nij,ni +1J | e { " Ar '^ l) | n i>j+1 , n m)i+1 >, (23) 

and Pi.j = 1 for all other choices (we used the notation At — ff)- 

Because both i and the j indices label occupation numbers, by using the Trotter-Suzuki 
(T-S) approximation we created a two-dimensional lattice from our real one- dimensional one. 
Every box of the original lattice is multiplied 2M times, and all the sets (j = 1,2M) 
are independent. Due to the original trace the j index must satisfy the periodical boundary 
conditions: n ii2 A/+i = 

We are able now to consider our quantum-statistical problem as a classical one, charac- 
terized by the 

Z^Y exp(-/3£ e M') (24) 

i n i,j} i d 
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partition function. In equation (|24f) we used the —fieij = InPij notation and every e^j is 
calculable from the four n it j, n i+ ij, Uij+i and n i+ ij + i occupation numbers. 

We have transformed our quantum-statistical problem in a classical one, considered in a 
space with increased dimensionality. The dimension indexed by % is called spatial, and the 
one indexed by j the Trotter or imaginary time direction. The presented method is totally 
equivalent with the Path Integral [6] formulation of the quantum-mechanical problem. (Each 
two lines in the Trotter direction corresponding to a time interval At in the Path Integral 
formalism.) 

The Monte Carlo simulation will follow now the known Metropolis algorithm. In the 
transformed, fl2H]) problem we have four-site interactions between the neighboring sites riij, 
rii+ij, Titj+i and rii+ij+i, with the values of i and j both even or odd (otherwise = 0). 
In the % and j space this interaction can be represented by a check-board pattern (Fig. 1), 
with interactions around the dark plaquettes. 

In many cases the P it j probabilities fl23|) are zero, and so the corresponding state is 
realized with zero probability. The condition for P it j ^ 0, is: 

n id + n i+1 ,j = n iJ+1 + n i+ld+1 . (25) 

One will realize immediately that in MC simulations the changes leading to acceptable 
configurations are quite limited. To save precious computer time and calculations, we must 
know from the beginning which changes will give nonzero transition probabilities. The 
acceptable changes (leading to nonzero transition probabilities ) will be: 

• the occupation numbers from the left side of a white plaquette are increased by unity, 
and the values from the right side decreased by unity 

• the occupation numbers from the left side of a white plaquette are decreased by unity, 
and the values from the right site increased by unity 

(These changes can be done only if the occupation numbers will not become negative, or for 
fermions if they are not bigger than one.) 

It is obvious that by satisfying initially the (^3) conditions for the whole system, and 
by considering only the above mentioned changes in the occupation numbers, the algorithm 
will lead to configurations keeping fl25|) . 

The Metropolis MC algorithm will be now as follows: 

• we consider an initial configuration so that condition ( p5|) holds for the whole lattice 

• we change the values of the occupation numbers around a white plaquette in the way 
described earlier 

• we calculate the change in the total energy of the system as AE = J^Aeij, the sum 
being done on the neighboring dark plaquettes 

• we accept this change with a probability 



• we continue the algorithm until thermodynamic equilibrium is approached 

• we collect periodically the relevant data. 

We define one MC step as 2ML trials of changing the configuration of the system. 



III. DETERMINING THE RELEVANT PHYSICAL QUANTITIES 

The relevant physical quantities are determined after the equilibrium dynamic is ap- 
proached, by averaging over many MC steps (usually of order 10 5 ). The E average energy 
and the C heat-capacity is computed by using the known equations: 

E = (27) 

ri _ 1 d 2 InZ /qq\ 

Keeping in mind that the factors are also functions of f3 from (|24]) we get: 

E=(j2e^ + ^E d -^ (29) 
C^(((ge, + ^^))-((ge, + /3g^V 



^(2E^ + /?E^) (30) 

■ h3 M 



In the above formula all the sums are over the dark plaquettes of the lattice and the averaging 
is an ensemble average, which in our QMC scheme is an average over different MC steps. 



IV. CALCULATION OF THE E^j FACTORS 

We will calculate the "energies" both for one particle and for many interacting 
fermions and bosons. For the later case we will assume that the interaction potential be- 
tween the particles in the same cell is two times stronger than the one between two particles 
in the nearest cells. We will also assume that there is no external potential energy, thus 
W t = 2t. 

A. One quantum particle 

The Hamiltonian of the particle depends on the kinetic energy only: 

H = ~t E(c+c, +1 + cf +1 c t ) + 2t ■ I> (31) 

i=l i=l 

We use the symmetric form: 
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Hi — —t • (c+q+i + cf +1 Ci) + t ■ hi + t ■ h i+ i where i — 1, L (32) 

L L 

if = £ifi = -t£di (33) 

i i=l 

6i = {cjc i+l + c+ +1 Ci) - hi- n i+1 , where i — l,L. (34) 
One can write the values as: 

I ( °° / (3t\ k 1 - \ 

e ij = ~p ■ ln I E (^J ' < > j (35) 

eij depends only on the occupation numbers riij, rij+ij, riij+i, n i+ ij + i, thus 
For one particle 

L 

E n M = 1 ' ( 36 ) 
i=i 

and only the terms with rijj + n i+ ij = n^+i + n i+1 J+1 can be different from zero: 
e^(0, 0, 0, 0), e^(l, 0, 1, 0), ^(0, 1, 0, l), 6^(1, 0,0, 1) and e^(0, 1, 1, 0). The e^(0, 0, 0, 0) term 
is zero, and by using the fact that the Oi operators are hermitian it is easy to prove: 

e iJ (l,0,l,0) = e i , J -(0, 1,0,1) 

e M (l,0,0,l) = e M (0, 1,1,0). (37) 

Using the immediate 

Oi|l,0 > = 1 • |0,1 > -1 • 1 1, > 

Oi|0,l > = 1 • |1,0 > -1 • |0,1 >, (38) 
equations, the two nontrivial terms can be calculated by a recursion formula: 

< l,0|6f|l,0 > = < 1, 0|6f _1 |0, 1 > - < l,0|O? _1 |l,0 > 

< l,0|Of|0,l > = < l,0|Of- 1 |l,0> - < l,0|Of _1 |0,l > . (39) 

Let us denote by a k =< l,0|Of|l,0 > and b k =< l,0|Of|0,l >, and rewrite the above 
recursion: 

Ofc = bh-i — O/t-i 

b k = a k -i - b k -i. (40) 

Because < 1, 0|1, >= 1 and < 1, 0|0, 1 >= the first terms of these series are a = 1 and 
b = 0. For k > we get the analytical forms: 

a k = (-l) k -2 k - 1 

b k = (-l) k+1 ■ 2 k -\ (41) 
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Substituting these into the expression of the terms ( j35|) we get: 

t; I vt — m | cosh 
(3\M \ \M 



MW,0) = ^(§-ln(cosh(^V (42) 



e,,(l,0,0,l)4g-ln(s 1 nh(§))). (43) 



B. Many interacting fermions 



The Hamiltonian in this case is: 

L 



H = ~tJ2 (cjfic* + cp j+1 ) + 2t ]T cp 3 + ^iE • ( 44 ) 

j'=i j'=i i=i 

Let us assume that the particles have charge q and there are rejective forces between 
them. The interaction is only between the nearest cells, because spin-less fermions can 
not be simultaneously in the same cell. (This is the manner we impose the antisymmetric 
wave-function for fermions). The interaction potential is V\ 

Vi — j— • — i (45) 

where a is the width of the cells. 
The Hamiltonian terms are 

Hi = -t(cfci +1 + cf +1 di - hi- h i+ i - 2bhi ■ h i+1 ), (46) 

with b = -y- The occupation numbers for spin- less fermions can be either zero or one, thus 
eij has three different values: (1, 0, 1, 0), e^l, 0, 0, 1) and ey(l, 1, 1, 1). The first two were 
already calculated. 

For the calculation of ey(l, 1, 1, 1) we show that 

6?|l,l>=-2(l + 6)6f- 1 |l,l>, (47) 

and get: 

< 1,1|6*|1,1 >= [-2(1 + 6)]*. (48) 
Substituting this in the expression of the energy: 

6^(1,1,1,1)= \ } . (49) 
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C. Many interacting bosons 



In this case more particles can be in the same cell, so the interaction between particles 
in the same cell has to be considered too. The Hamiltonian has the following form: 

L L L L 

# = -* E + cfc j+1 ) +2tJ2 %Cj + V 1 J2 n 3 h 3+1 + V -J2 " !)■ (50) 

j=i j=i j=i j=i 

Let us assume V = 2V±. In our standard notations: 

Hi = -t[cjc i+1 + cj +1 Ci - hi- h i+1 - 2bhih i+1 - bh^hi - 1) - bh i+1 (h i+1 - 1)]. (51) 

In case of N bosons we have two conditions: 

L 

E n M = N > 
i=l 

n id + n i+ ij = n id+1 + n i+1 j +1 . (52) 

Due to the % «-> % + 1 and j j + 1 invariance of the values we consider only the 
dark squares with the biggest occupation number in the left bottom corner. For simplicity 
we introduce the following notation 

< nij, ni + ij\Oi\nij + i, >= (riij, n i+ ij, ra i+lj _,- +1 ) fe , (53) 

and calculate the recurrence formula for a given n = n,ij+n i+ ij = nij + i+n i+ ij + i occupation 
numbers (n £ [0, N]). If n is odd and n > 1: 

fn + 2l + l n-2l- \ n + 2k + 1 n - 2k - l\ 
{ 2 ' 2 ' 2 ' 2 ) k ~ 



(n + 2k + 3)(n-2k- 1) /n + 2/ + 3 n - 2/ - 3 n + 2A; + 1 n-2jfe-l\ 



fc-i 



(n + 2A; + l)(n-2A; + l) f n + 21 + 1 n - 21 - 1 n + 2k + 1 71-2/5-1' 



2 A_x 



m u ,,Jn + 2l + l n-2l-l n + 2k + l n-2k-l\ 
n(l + Kn - 1)) 5— , — , j ' 2 J k " (54) 



The possible values of / and k are / G 0, ! ^ and fee [—/ — !,/] 
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If n is even 



n n n n 
2 ' 2 2 ' 2 /k 



In + 2k + 2n - 2k (n + 21 n - 21 n + 2k + 2 n - 2fc - 2\ 



2 V 2 ' 2 ' 2 



ln + 2kn-2k + 2 f n + 2l n - 21 n + 2k - 2 n - 2k + 2 



fc-i 



ra(l + &(n- 1)) 



'n + 2/ n-2Z n + 2A; n-2k s 



(55) 



where / G 



- 



and fe G [— /, /]. 

This equations carries all the necessary information for calculating numerically the 
energies. 



V. TEST FOR ONE PARTICLE IN AN INFINITE DEEP POTENTIAL VALLEY 

We verified our QMC algorithm by considering first one quantum particle (an electron) 
in an infinite deep potential valley and in contact with a heat-bath at temperature T. This 
simple quantum-statistical problem is easily computable. The energy levels are given by the 
well-known formula: 

B "=8^' < n=1 - 2 >• (56) 

(ah is the width of the valley, and was taken 3 • lCT 9 m) The expectation value of the energy 
is numerically computable as: 

< E >f ~ • (57) 

One can also easily compute theoretically the heat-capacity, of the system by numerically 
derivating < E > t as a function of T. 

Considering a series of simulations with L — 12 and p = (3t/M = 0.1 we got the 
< E >qmc values plotted with dots in Fig. 2a. This result is presented in comparison 
with the < E > t (T) theoretical curve. On Fig. 2b we plotted together the theoretically 
computed heat-capacity (dashed line), the one obtained by QMC applying formula (|30|) 
(filled circles) and the one obtained from the 

C = A< f r >0MC . (58) 
formula (empty circles). The correspondence is satisfactory. 
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VI. PROBLEM WITH ONE QUANTUM PARTICLE ON A RING 



After the promising results for the infinite deep potential valley it is quite frustrating 
to realize that in this case we cannot solve the original quantum-statistical problem. In 
the QMC method imposing periodic boundary conditions in the spatial direction will define 
a problem for a quasi-free particle in constant potential field. The reason for this is that 
in the QMC formalism we loose the wave-function and remain only with the occupation 
probabilities. In this manner it is impossible to impose the necessary closing conditions for 
the wave-function and it's derivative. From the symmetry of the problem we get only the 
condition that the occupation numbers have translational invariance. This defines a problem 
for a free quantum particle (i.e infinite wide potential valley). We expect in this manner to 
find 



< E > t = f , (59) 
C = f , (60) 



which are the known results for a free particle in contact with a heat-bath at temperature 
T. In contrast, by imposing the ring-like geometry in the quantum- mechanical problem, the 
closing condition for the wave-function would yield the 



h 2 n 2 

£ "=WI^ < B = 1 ' 2 ' 3 »• < 61) 

energy levels and a different < E >' t (T) dependence. 

The previous predictions are totally supported by our simulation results. We considered 
a quantum particle with the mass of an electron, aL = 3 ■ 10~ 9 m and L = 18. On Fig. 3a we 
plotted the < E >qmc simulation points together with the expected < E >= kT/2 curve 
(continuous line) and the solution with the discrete (|6T| ) energy levels for the quantum states 
(dashed line). The simulation results are convincingly supporting the free quantum particle 
solution. The heat-capacity obtained from formula fl30l) fluctuates around the expected k/2 
value (Fig. 3b). 

Through this simple exercise one can immediately realize a very important limitation 
for the QMC methods. By loosing the *Sf wave-function, closing or boundary conditions for 
\I/ and it's derivative are replaced with continuity conditions for | ^ | 2 and particle flux j. 
The two problems however are not equivalent, and important differences can be obtained 
for several problems. 



VII. MINUS-SIGN PROBLEM FOR INTERACTING BOSONS 

For bosons where the occupation numbers inside one cell can be arbitrary high an im- 
mediate difficulty arises when calculating the (|23| ) P^ transition probabilities. For large 
occupation numbers the P^ values become negative, and thus the QMC method is inap- 
plicable. This problem is similar with the well-known minus-sign problem for the Hubbard 
model H . The origin of this non-physical situation is the used Trotter-Suzuki approxima- 
tion and thus the finite value of p = (3t/M. As the value of p is lower more bosons can be 
considered inside one spatial cell without getting minus sign for P^. On Fig. 4 as a function 
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of p we plotted the maximal boson number per cell for which the Pij factors are all positive. 
(We considered here free bosons with the mass and charge of an electron interacting with 
a repulsive Coulomb potential, and a = 1.66 • 10~ 10 m.) From this figure we learn that in 
the limit p — > the minus-sign problem disappear. Unfortunately in this limit the M value 
(lattice size in the imaginary Trotter direction) has to be infinite. The practical solution 
is of course by fixing the maximal number of bosons in the simulation. This leads to a 
maximal value of p, and thus a finite M value. In this manner the minus-sign problem will 
be eliminated. The technical problem we will face now is that even for small number of 
bosons the required M value is large and this limits seriously the applicability of the QMC 
method for the proposed problem. 



VIII. UNPHYSICAL FRUSTRATION FOR INTERACTING FERMIONS 

In MC-type simulations the usual method for calculating the heat-capacity is by the use 



of the fluctuation-dissipation theorem. In our case this leads to equation fl30D , and we get 
a C\ value. An alternative method for calculating the heat-capacity (C2) could be by the 
definition of this quantity, i.e. by using equation (|58|). 

Considering the problem of interacting fermions (spin-less electrons) on a ring-like ge- 
ometry we proposed to determine the specific heat of the system both as a function of 
temperature and particle number in the system. We considered repulsive Coulomb inter- 
actions, a = 1.66 ■ 10 _10 m, and L — 18 . In order to get confidence in our calculations we 
used both of the methods. The obtained results were in agreement for low and high fermion 
number, and totally contradictory for medium occupation density. The logarithm of the ra- 
tio (Ci/C 2 ) of the two specific heat is plotted in Fig. 5. Verifying and testing the result in 
several aspects we got confidence in the obtained differences. The explanation of this strange 
behavior lies in the space discretization and the way of handling the exclusion principle in 
QMC simulations. The combined effect of this two introduces a non-physical frustration in 
the system. This frustration is maximal for half-filled cells and leads to results similar with 
the ones obtained in spin glasses below the freezing temperature. In MC simulations for spin 
glasses the difference between the specific heat determined in the two different manner is 
known ||, and indicates that thermal equilibrium is not completely reached. We encounter 
a similar situation also here. 

To better understand our previous arguments let us follow what is happening in our 2D 
lattice following the particles paths. Each fermion is represented by a continuous series of 1 
values in the Trotter direction which defines the "path" (or the "world-line") of the particle. 
For fermions these paths cannot intersect each other and can change maximally one cell size 
for 1 step in the Trotter dimension. A possible path configuration for two fermions is shown 
as an example in Fig. 6. If we consider more than two fermions, due to the introduced space 
discretization and the restrictive condition for fermion numbers inside one cell many paths 
become frustrated. A typical situation for three fermions is shown in Fig. 7. In this setup 
the path for the fermion in the middle is frustrated in our QMC algorithm. By a no energy 
change this path can fluctuate between the two neighbors. The energetically stable position 
is in the middle, but due to our space discretization there is no cell there. We expect that 
this effect is maximal around the half-filling and disappears for less than three fermions or 
for complete filling. This prediction is in agreement with the observed results. 
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The solution for eliminating this unphysical frustration is by considering a fine space 
discretization where the frustration is minimized. Unfortunately this will increase again the 
lattice size and create technical problems. 

IX. CONCLUSIONS 

Considering the quantum statistical problem of interacting particles on a ring-like geom- 
etry and in an infinite deep potential valley we realized several limits of the QMC method. 

• We loose the concept of wave-function and the method becomes inappropriate for ring- 
like geometries. Imposing periodic boundary conditions will lead to the problem of 
particles in a spatially constant potential. 

• For interacting bosons we found that due to the considered Trotter-Suzuki transfor- 
mation a minus-sign problem appears. As the Trotter dimension is increased (the 
Trotter-Suzuki approximation is improved) more and more bosons per spatial cell can 
be considered without getting into the minus-sign problem. This problem can be 
eliminated by increasing the lattice size, which leads to increased computational time. 

• For interacting fermions the combined effect of spatial discretization and the applica- 
tion of the exclusion principle creates an unphysical frustration in the system. In this 
way thermal equilibrium is hard to reach and problems similar to the one obtained in 
MC methods for spin-glasses are encountered. The effect of this unphysical frustration 
can be minimized by further increasing the spatial cell numbers, which creates again 
technical difficulties. 

In conclusion, for a successful application of the QMC method on the proposed problem 
one has to keep in mind and eliminate the above limitations. 
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FIGURES 
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FIG. 1. The characteristic chess-board lattice for the QMC simulations 
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FIG. 2. QMC results in comparison with the exact theoretical results for an electron in a ID 
potential valley with width aL = 3 • 10 -9 m. Fig. 2a presents the energy values of the system as a 
function of the heat-bath temperature. The dashed curve is theoretical, the filled circles are QMC 
results. Fig. 2b presents the heat-capacity of the system as a function of the temperature. The 
dashed curve is the exact theoretical result, filled circles are QMC results obtained by (|30|), and 
the empty circles are results obtained from (j58|). (L = 12 and p = 0.1) 
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FIG. 3. QMC results in comparison with theoretical results for an electron on a ring (length 
of the ring: 3 • 10 _9 m). Fig. 3a show the energy values as a function of the heat-bath temper- 
ature. Filled circles are QMC simulation results, dashed curve is the desired result for the right 
quantum-mechanical problem, and the continuous curve presents the theoretical result for a free 
electron (i.e. infinite large ring). On Fig. 3b we present the heat-capacity of the system as a 



function of the temperature. Circles are QMC results, obtained by (30), and the line represent the 
theoretical k/2 value for a free particle. (L = 18 and p = 0.2) 
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FIG. 4. Maximal bosons number per cell (n), for which the minus-sign problem does not appear 
as a function of the p parameter. We considered bosons with the mass and charge of the electron, 
interacting with repulsive Coulomb potential, (a = 1.66 • 10~ 10 m and L = 18) 
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FIG. 5. Logarithm of the ratio for the heat-capacity calculated by two different ways {C\ and 
C2), as a function of particle number per cell. (Interacting fermions with mass and charge of an 
electron, Coulomb potential, aL = 3 • 10 _9 m and L = 18.) 




FIG. 6. A specific path configuration for two fermions 
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FIG. 7. Example of frustration for three fermions 
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